Normal-state magnetic susceptibility in a bilayer cuprate 
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The magnetic susceptibility of high- Tc superconductors is investigated in the normal state using a 
coupled bilayer model. While this model describes in a natural way the normal-state pseudogaps 
seen in c-axis optical conductivity on underdoped samples, it predicts a weakly increasing suscepti- 
bility with decreasing temperature and cannot explain the magnetic pseudogaps exhibited in NMR 
measurements. Our result, together with some experimental evidence suggest that the mechanism 
governing the c-axis optical pseudogap is different from that for the a-b plane magnetic pseudogap. 
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Normal-state pseudogaps of underdoped samples have 
been one of the important subjects in the study of high- 
Tc cuprates in the past few years (l). Experimentally 
strong decrease of NMR Knight shift and spin-lattice re- 
laxation rate with decreasing temperature has been ob- 
served in various normal-state underdoped samples below 
some cross over temperature T* [^j-^. These decreases 
have often been considered in terms of the opening of a 
"spin gap" due to the antiferromagnetic (AF) spin fluc- 
tuations or, alternatively, the notion of "preformed 
Cooper pairs" in the case of significant order parame- 
ter phase fluctuations [|lO|. Measurements on dc resis- 
tivity for underdoped samples also reveal that the in- 
plane resistivity pab{T) deviates from a linear in tem- 
perature behavior at some temperature consistent with 
the crossover temperatures T* in NMR experiments [0 . 
On the other hand, the c-axis optical conductivity cTc{ijj) 
|p"2|-p^ exhibits a striking gap-like depression at low fre- 
quency and the c-axis resistivity Pc(T) displays a upturn 
semiconducting feature |pT[| in the normal-state under- 
doped samples, characteristic of the formation of some 
kinds of gap. 

A central question in the study of pseudogaps is - are 
the pseudogaps seen in the c-axis optical conductivity 
related at all to the magnetic pseudogaps seen in NMR 
experiments? Or, could there be two different mecha- 
nisms? While there is currently no consensus on this 
point, we begin by comparing first some of the properties 
of these two pseudogaps seen in experiments. Consider 
the magnitude and anisotropy of these two pseudogaps: 
(i) Magnitude ~ Although one does not measure the mag- 
netic pseudogap directly, an analysis done by Williams et 
al. of Knight-shift data based on quasiparticle exci- 
tation spectrum given by i?k = [e^ -I- A^]^/^ where Ak 
is the normal-state pseudogap shows that the magnitude 
of magnetic pseudogap is strongly dependent on doping 
or carrier density (larger gap with lower doping) and 
scales with the crossover temperature T* . In contrast, 
the c-axis optical conductivity measured by Homes et al. 
@ on YB2C3O6+X (YBCO) at different doping x in- 
dicates that the frequency ranges where conductivity is 
suppressed (direct measurement for the size of the gap) 
are almost independent of doping. Moreover, a compari- 



son § of YBCO and Pb2Sr2(Y,Ca)Cu308+5 (PSYCCO) 
data gives that the optical pseudogap does not scale with 
T*. (ii) Anisotropy - NMR experiments reveal that 
the magnetic pseudogap has d-wave like anisotropy (con- 
sistent with the angular resolved photoemission spec- 
troscopy (ARPES) ijr^ result), while the flat spectral 
weight seen within the pseudogap region in c-axis op- 
tical data implies that the optical pseudogap is isotropic 
(s-wave like). The above comparison seems to suggest 
that the magnetic pseudogap is intimately connected to 
the superconducting gap and reflects the pairing mech- 
anism, while the optical pseudogap is due to a different 
mechanism. 

Recently Atkinson and Carbotte and with Wu 
p9| , p0[ have used a coupled bilayer model to calculate 
various c-axis properties of the high Tc cuprates. Within 
this model, a minimum energy difference A between two 
non-degenerate bands (associated with two layers cou- 
pled by a single-particle hopping tj_) is introduced. This 
band gap A - turns out to be the pseudogap seen in 
the c-axis optical conductivity - and gives the suppres- 
sion of the joint density of states in the interband tran- 
sitions. The key feature of this coupled bilayer model is 
that the c-axis optical conductivity can be generally sep- 
arated into intraband (proportional to f^) and interband 
(proportional to t'^) parts. Therefore, when t± is small 
which is the condition appropriate to the underdoped 
regime pl[ , interband contribution dominates over the 
usual intraband contribution and, as a result, the pseu- 
dogap is visible (non-Drude like behavior). In contrast 
in the larger t± optimally-doped or overdoped regime, 
intraband dominates and one retains the more usual no- 
pseudogap Drude-like behavior. 

It is shown in Refs. [^,^ that the coupled bilayer 
model explains in a natural way the pseudogap seen in 
the c-axis optical conductivity for underdoped, and gives 
an explanation of why it is not seen in optimally-doped 
or overdoped samples (see Fig. ^). Both temperature and 
frequency dependence of the entire spectrum seen in ex- 
periments |p^ , p^ are well described by this model. It is 
thus not surprisingly that the semiconductor-like upturn 
feature at low temperature seen in the c-axis dc resis- 
tivity is also understood in this model EG]. Since this 
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FIG. 1. Normal-state c-axis conductivity as a function of 
photon energy for a plane-chain bilayer at different temper- 
atures with a small t± = 2meV. The minimum energy dif- 
ference between the two bands for a given k in the first 2D 
Brillouin zone is A = 20meV (taken from Fig. 1 in Ref. @). 



model emphasizes the band structure correlation arising 
from two non-degenerate coupled layers, a simple test can 
be immediately done for the single-layer La2-2:Sr2:Cu04 
(LSCO) which is found by Basov et al. to exhibit no 
pseudogap in the c-axis optical conductivity or a very 
weak pseudogap with large energy scale O.leV by 
Uchida et al. [ p3[ . The latter case accompanied with 
a semiconducting upturn feature of Pc{T) strongly 
suggests that while there is only one conducting layer, 
in practice there are two bands separated by a large en- 
ergy difference A ^ O.leV in LSCO. We note that in this 
case Startseva et al. [ p5| find a strong in-plane pseudo- 
gap with smaller energy scale indicating no correlation 
between in-plane and out-of-plane energy scale for the 
pseudogap. 

In the following, we use the bilayer model to calculate 
the magnetic susceptibility. We consider the spin raising 
(+) and lowering (— ) operators for layer i {i—1,2) defined 



by 



and 



k ''i,k-|-q,i''^k,T- 



The transverse (parallel to a-b plane) magnetic suscepti- 
bility Xij which we shall concentrate on, is then given by 

pel 



dre-"-(T.ar^(r)4_q(0)), (1) 



where /ig is the Bohr magneton, Xii,X22 denote the 
in-plane spin correlation, and Xi2,X2i denote the out- 
of-planc (or cross) spin correlation. Due to its local 
nature, an NMR experiment is able to study the spin 
susceptibility Xij (Qi ^) individually. For example, the 
Knight shift {Ks) for an atom in layer i is directly 
proportional to the real part of the planar spin sus- 
ceptibility, Ks{T) ~ Xii(q = OjW = 0) and the spin- 
lattice relaxation rate (per temperature unit) is related 



to imaginary-part susceptibility over a weighted sum, 
(l/TiT), - EqJ^(q)l'xr,(q:'^ - 0)/u;. While, the in- 
formation contained in the cross spin susceptibility Xi2 
can be probed by a spin echo double resonance (SEDOR) 
experiment (see, for example, Ref. p^). 

In terms of a free electron system, the susceptibility 
(111) can be reduced to 



Tr[G'o(k, w„)7jGo(k -f- q, w„ - iw„)7j 



tracing over the two bands, where 71 = 
are the vertices and 
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(2) 
and 

(3) 



is the Green's function matrix. The parameter ^.^ is the 
band structure for isolated layer i and t is the coupling 
between the two layers. Considering only the q — > case, 
the trace operator in (|^) enables one to calculate Xij iii 
a convenient frame in which the Green's function matrix 
is diagonal. Consequently, 



Xij (q = 0, iw„) = -2^ll XI ^ XI 



Tr[G'o(k, it'„)7jGo(k, if„ - ia;„)7j 



(4) 



where 



Gg ^(k, w„) = 

ivn - e+ + ir+sgn(i/„) 

iVn 





£- + ir_sgn(t'„) 



(5) 



with T± introduced as the total scattering rate in band 
±. The two renormalized bands are 



and the rotated vertices are 
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where an = (^1 - e_)/ (e+ - e_), a22 = (e+ - 6)/(e-i- - 
e_), and ai2 = |t|/(e+ — £-)• Substituting (||) and (|^) 
into (^), we obtain 

Xn(«t^n) = -2AieX;^X [("ll'^ll^il +"22G22G22) 



X22(jt^n) 



^12(^11^22 + G'22G'ii 

-2/^' E ^ E [("22GnG;i + a?iG22G'22) 
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>^12(G'iiG22 + G22G'ii 



(8) 



where Gu = Gji(k, ii/„) and G-j = Gii(k, w„ - iujn) are 
the elements of (diagonal) Green's function matrix in (|^). 
Clearly for each Xij in l^)? the first term corresponds 
to an intraband contribution and the second term cor- 
responds to an interband contribution. One sees that 
the in-plane xu and cross Xij (* 7^ j) have equal but 
opposite-sign interband contributions. 

The frequency sum in (H) can easily be transformed 
into an integral such as (using (w) 



da; 



2r 



_^27r [(x-e,)2 + r2][(:r-c.-e,)2+r2] 



/(x)(a; - cj - ej) + f{x - uj){x - e^) 



(9) 



where /(a;) is the Fermi distribution function and we 
have redefined ei = e+ and €2 = The scattering 
rates are simply assumed to be identical for both bands, 
r_|_ = r_ = r. To obtain analytical results, we consider 
the two isolated layers to have the following simple band 
structures 



2m 



£.2{kx,ky) = — fcjj - 



(10) 



ky)^^^ is the 2D layer momentum, fj, is 



where fc|| = (fc^ 

the chemical potential, and A corresponds to the energy 
difference between these two layer bands (in general, A 
corresponds to a minimum energy difference between two 
bands). For YBCO, we consider the Cu02 plane as layer 
1 which has generally more (hole) carriers (or less elec- 
trons) than the CuO chain which is layer 2. While the 
band structure (^) is oversimplified, it has captured the 
essential physics built in this coupled bilayer model. Tak- 
ing into account the periodicity of crystals, the layer-layer 
couphng is chosen to be t{kz) = 2t± cos{kzd/2), with the 
parameters t± the coupling strength and d the spacing 
between two layers. When t± is small {\t\ ^ A), one may 
expand an ~ 1 - t^/A^, 022 ^ f/A'^, and ai2 ~ \t\/A. 
This in turn implies that for the planar spin susceptibil- 
ities xii and X22, the intraband term is of order 1 and 
the interband term is proportional to i^, while X12 has 
both intraband and interband terms proportional to t\. 

With the above simplifications, the real part of the 
magnetic susceptibilities in (^) at zero frequency can be 
found to be (keeping only second order terms in t±) 
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FIG. 2. Temperature dependence of the theoretical planar 
magnetic susceptibility Xn(0) given in (|ll|). The curves from 
top to bottom correspond to fi = 400, 200, and lOOmeV with 
fixed A = 20meV, t± = 3meV, and r(r) = 20meV at lOOif 
and is linear in temperature. 
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where iV(0) = m/Trh^ is the 2D (constant) density of 
states for free electron gas and we have defined 



arctan 
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arctan 



Ai - A 



7(T) ^In 



(m- A)2 + r2 
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In deriving (|T^), we have made use of the fact that 
—df{x)/dx sa 6{x), which is appropriate for not too 
high temperatures. In Fig. ^ we plot the temperature- 
dependent layer-1 Knight shift (proportional to xii(O) 
in (pd])) using the parameters which is appropriate to 
and has successfully described the c-axis optical pseudo- 
gap for underdoped YBa2Cu306.6: A* = 100 to 400meV, 
A = 20meV, t± = 3meV, and T{T) = 20meV at lOOi^ 
and is linear in temperature. It is clearly shown in Fig. ^ 
that the Knight shift increases only slightly as temper- 
ature decreases and docs not exhibit the decreasing fea- 
ture characterizing the magnetic pseudogap seen in ex- 
periments. Since in YBCO, A <^ fi, X22(0) is found to be 
similar to xii(O) of Fig. |2| Apart from the magnitude of 
tj_, it can be shown quite generally that when T(T) < /z, 
/3 - 7 -r^/[/i^(Ai^ -I- r^)] which is small. This means 
that the last term in Xii(O) o^' X22(0) in (|ll|) does not 
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contribute much and one expects a similar temperature- 
dependent Knight shift as given in Fig. ^ even for over- 
doped YBa2Cu307 for which t± is much larger. Naively 
if one applies the result x'n (0) in (|Tl| ) to both YBCO with 
smaller A = 20meV and LSCO with larger A ^ O.leV 
and assumes /i and F identical for both cases, one expects 
more significant increasing Knight shifts with decreasing 
temperature in overdoped LSCO compared to those in 
overdoped YBCO which seems consistent with experi- 
ments Hi. 

In contrast to the c-axis conductivity where the inter- 
band term (cx t'^) dominates over the intraband term 
(oc i^) in the small t± underdoped case, for the planar 
spin susceptibility, the interband term (oc t^) is almost 
canceled by an equal contribution arising from intraband 
terms. Consequently, the planar spin susceptibility is not 
very much dependent on the interband transition (or, 
equivalently, the magnitude of tj_), and, as a result, also 
not very dependent on temperature. Thus the Knight 
shift derived from the simple coupled bilayer model has 
only a weak temperature dependence mainly due to im- 
purity scattering effect (i.e., F(T)). 

In principle one can measure directly the effect of t± 
with SEDOR . An effective coupling a exists between 
spins on the two layers ali •l2- In terms of our model a ~ 
hA^Xi^/ ^J''B where A is the hyperfine coupling constant. 
In order to obtain an estimate for a we take A to be the 
isotropic transferred hyperfine constant for Cu atoms in 
the CuOz planes of YBCO where A/^^j = 82 kOe and 
take iV(0) = 2 states/eV-Cu. Evaluating Xi2 at T = 300 
K with II = 100 meV yields Xu/l^B^i^) = 0.0016. Thus 
the effective coupling rate is a ~ 15 s~^ which is about 3 
orders of magnitude too small to observe experimentally 

While the theory given in Eq. ( pi] ) does not describe 
the magnetic pseudogaps for underdoped samples, sev- 
eral predictions and conclusions can be drawn in con- 
nection with experimental data. First, there exists a 
mechanism (not the interband effect) which operates 
most profoundly in the underdoped samples and leads 
to the suppression of Knight shifts (magnetic pseudo- 
gaps) . Secondly, this mechanism probably contributes to 
the pseudogap effects seen in all the a-b plane measure- 
ments, including the charge-associated in-plane dc resis- 
tivity pab{T), ARPES, in-plane scattering rate 1/t{u!,T) 
in the infrared reflectivity, and the spin-associated NMR 
measurements and is of only secondary importance in c- 
axis transport. Thirdly, the interband effect we propose 
should dominate in the c-axis transport and is strongly 
associated with the pseudogap effects seen in c-axis opti- 
cal and dc resistivity measurements. Fourth, as indicated 
in experiment on YBa2Cu408 [Q, the magnetic pseudo- 
gap mechanism has a much stronger effect on the Cu02 
plane (deeper normal-state Knight-shift drop) than CuO 
chain (weaker or almost no drop). 

In summary, while at present one cannot rule out that a 
single mechanism may govern both the c-axis optical and 
a-b plane magnetic pseudogaps, we have shown through 



our theoretical studies as well as with experimental data 
that a two-mechanism picture is favored. 

The authors thank Bozidar Mitrovic, Tom Timusk, 
and Tatiana Startseva for very stimulating discussion. 
This work was supported by Natural Sciences and Engi- 
neering Research Council (NSERC) of Canada and Cana- 
dian Institute for Advanced Research (CIAR). 



[1] See, for example, Barbara Goss Levi, Physics Today 

49(6), 17 (1996). 
[2] M. Takigawa, A.P. Reyes, P.C. Hammel, J.D. Thompson, 

R.H. Heffner, Z. Fisk, and C. Ott, Phys. Rev. B 43, 247 

(1991) . 

[3] R.E. Walstedt, R.F. Bell, L.F. Schneemeyer, J.V. 
Waszczak, and G.P. Espinosa, Phys. Rev. B 45, 8074 

(1992) . 

[4] M. Bankay, M. Mali, J. Roos, and D. Brinkmann, Phys. 

Rev. B 50, 6416 (1994). 
[5] R. Stern, M. Mali, I. Mangelschots, J. Roos, D. 

Brinkmann, J.-Y. Genoud, T. Graf, and J. Muller, Phys. 

Rev. B 50, 426 (1994). 
[6] Ya-Wei Hsueh, B. W. Statt, M. Reedyk, J. S. Xue, and 

J. E. Greedan, Phys. Rev. B, submitted. 
[7] A. J. Millis, H. Monien, and D. Pines, Phys. Rev. B 42, 

167 (1990). 

[8] H. Monien, D. Pines, and M. Takigawa, Phys. Rev. B 43, 
258 (1991). 

[9] A.V. Chubukov, D. Pines, and B.P. Stojkovic, preprint. 
[10] V. J. Emery and S. A. Kivelson, Nature 374, 434 (1995). 
[11] K. Takenaka, K. Mizuhashi, H. Takagi, and S. Uchida, 

Phys. Rev. B 50, 6534 (1994). 
[12] C. C. Homes, T. Timusk, R. Liang, D. A. Bonn, and W. 

N. Hardy, Phys. Rev. Lett. 71, 1645 (1993). 
[13] C. C. Homes, T. Timusk, D. A. Bonn, R. Liang, and W. 

N. Hardy, Physica C 254, 265 (1995). 
[14] D. N. Basov, R. Liang, B. Dabrowski, D.A. Bonn, W.N. 

Hardy, and T. Timusk, Phys. Rev. Lett. 77, 4090 (1996). 
[15] M. Reedyk, T. Timusk, Y.-W. Hsueh, B.W. Statt, J. S. 

Xue, and J. E. Greedan, preprint. 
[16] G.V.M. WiUiams, J.L. Tallon, E.M. Haines, R. Michalak, 

and R. Dupree, Phys. Rev. Lett. 78, 721 (1997). 
[17] See, for example, A. G. Loeser et al, Science 273, 325 

(1996). 

[18] W. A. Atkinson and J. P. Carbotte, Phys. Rev. B 55, 
3230 (1997). 

[19] W. A. Atkinson, W. C. Wu, and J. P. Carbotte, Phys. 

Rev. Lett., submitted. 
[20] W. C. Wu, W. A. Atkinson, and J. P. Carbotte, Phys. 

Rev. B, submitted. 
[21] Y. Zha, S. L. Cooper, and D. Pines, Phys. Rev. B 53, 

8253 (1996). 

[22] D. N. Basov, H.A. Mook, B. Dabrowski, and T. Timusk, 

Phys. Rev. B 52, R13141 (1995). 
[23] S. Uchida, K. Tamasaku, and S. Tajima, Phys. Rev. B 

53, 14558 (1996). 



4 



[24] G. Boebinger et al, Phys. Rev. Lett. 77, 5417 (1996). 

[25] T. Startseva et al, preprint. 

[26] Sec, for example, S. Doniach and E.H. Sondheimer, 
Green's Functions for Solid State Physicists (W.A. Ben- 
jamin, Reading, Massachusetts, 1974), p. 162. 

[27] C. Slichter, Principles of Magnetic Resonance, 3rd ed. 
(Springer- Verlag, New York, 1990). 

[28] T. Nakano, M. Oda, C. Manabe, N. Momono, Y. Miura, 
and M. Ido, Phys. Rev. B 49, 16000 (1994). 

[29] B.W. Statt, L.M. Song, and C.E. Bird, Phys. Rev. B, 
May 1997 issue. 



5 



